Unconventional critical behavior of fermions hybridized with bosons 



A.S. Alexandrov 

Department of Physics, Loughborough University, Loughborough LE11 3TU, United Kingdom 

A phase transition into the condensed state of fermions hybridized with immobile bosons is ex- 
amined beyond the ordinary mean-field approximation (MFA) in two and three dimensions. The 
hybridization interaction does not provide the Cooper pairing of fermions and the Bose-condensation 
in two-dimensions. In the three-dimensional (3D) boson-fermion model (BFM) an expansion in 
the strength of the order parameter near the transition yields no linear homogeneous term in 
the Ginzburg-Landau-Gor'kov equation. This indicates that previous mean-field discussions of the 
model are flawed in any dimensions. In particular, the conventional (MFA) upper critical field is 
zero in any-dimensional BFM 
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Soon after Anderson pj and Street and Mott 0] in- 
troduced localized electron pairs to explain some un- 
usual properties of chalcogenide glasses, a two component 
model of negative U centers coupled with the Fermi sea 
of itinerant fermions was employed to study supercon- 
ductivity in disordered metal-semiconductor alloys 0, Q . 
When the attractive potential U is large, the model is 
reduced to localized hard-core bosons spontaneously de- 
caying into itinerant electrons and vice versa, different 
from a non-converting mixture of mobile charged bosons 
and fermions 0, Later on the model was applied 
more generally to describe pairing electron processes with 
localization-delocalization Q, and a linear resistivity in 
the normal state of cuprates 8]. The model attracted 
more attention in connection with hi gh- temperature su- 
perconductors |3, CU EH EH BJl Ej> Ell- I n particular, 
Refs. (TUTU claimed that 2D BFM with immobile hard- 
core bosons is capable to reproduce some physical prop- 
erties and the phase diagram of cuprates. 

The model is defined by the Hamiltonian, 

H = E &cL^.<r + EoY, b l b <i+ (!) 

k,o-=T,i q 

giV _1/2 ^ (0k& q C_ k +q/2,t C k+q/2,i + H.C.) , 
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where £k = — 2t(cosfc a + cosfcj,) — fi is the 2D energy 
spectrum of fermions, Eq = As — 2/li is the bare bo- 
son energy with respect to the chemical potential fi, g 
is the magnitude of the anisotropic hybridization inter- 
action, 0k = </>-k is the anisotropy factor, and N is the 
number of cells. Ref. 0] argued that 'superconductiv- 
ity is induced in this model from the anisotropic charge 
exchange (hybridization) interaction (g<^k) between the 
conduction-band fermions and the immobile hard-core 
bosons', and 'the on-site Coulomb repulsion competes 
with this pairing' reducing the critical temperature T c 
less than by 25%. Also it has been argued [I^, that the 
calculated upper critical field of the model fits well the 
experimental results. 

This, as well as some other studies of BFM applied a 
mean-field approximation (MFA) to the condensed phase 
of BFM, replacing zero- momentum boson operators by c- 



numbcrs and neglecting the boson self-energy. MFA led 
to a conclusion that 'bosons exist only as virtual state' 
at sufficiently large boson energy Eq, so that BFM ex- 
hibits features compatible with BCS characteristics with 
a relatively small fluctuation region GiE2, and describes 
a crossover from the BCS-like to local pair behavior 
13 However, our study of BFM beyond 
MFA revealed a crucial effect of the boson self-energy 
on the normal state boson spectral function. The energy 
of zero-momentum bosons is renormalized down to zero 
at the critical temperature T = T c , no matter how weak 
the boson-fermion coupling and how large the bare boson 
energy are. As a result, the Cooper pairing of fermions 
via virtual unoccupied bosonic states is impossible, but 
it could occur only simultaneously with the Bose-Einstein 
condensation of real bosons in BFM. 

Here I show that there is no BCS-like condensed state 
in the two-dimensional model. More surprisingly, the 
mean-field approximation appears meaningless even in 
three-dimensional BFM because of the complete boson 
softening. The phase transition is never a BCS-like 
second-order phase transition. In particular, the con- 
ventional upper critical field is zero in 3D BFM. 

Replacing boson operators by c-numbers for q = in 
Eq.(l) one obtains a linearised BCS-like equation for the 
fcrmion order-parameter Ak, 



»2, 



EnN 



Ak' taah(g k >/2fc B r) 
2£ k < 



(2) 



with the coupling constant g 2 = g 2 (l — 2n B ), renor- 
malized by the hard-core effects. Using a two-particle 
fermion vertex part in the Cooper channel one can prove 
that this equation is perfectly correct even beyond the 
conventional non-crossing approximation [Tflj . Neverthe- 
less, the problem with MFA stems from an incorrect defi- 
nition of the bare boson energy with respect to the chem- 
ical potential, Eq(T). This energy is determined by the 
atomic density of bosons (n B ) as (Eq.(9) in Ref. gj) 

tanh-^_ = l-2n B . (3) 
2k B T K ' 

While Eq.(2) is correct, Eq.(3) is incorrect because the 
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boson self-energy £&(q, fl n ) due to the same hybridiza- 
tion interaction is missing. At first sight EH the self- 
energy is small in comparison to the kinetic energy of 
fermions, if g is small. However £;,(0, 0) diverges loga- 
rithmically at zero temperature [l3l ] , no matter how week 
the interaction is. Therefore it should be kept in the 
density sum-rule, Eq.(3). Introducing the boson Green's 
function 



D(q, fi„) 



1 - 2n L 



iQ n — Eo — £&(q, O n ) 
one must replace incorrect Eq.(3) by 



(4) 



N 



n B , 



(5) 



where r = +0, and Q n = 2-nk B Tn (n = 0,±1,±2.. 
The divergent (cooperon) contribution to Sb(q,O n ) is 
given by 



2i(q,nn) = -^E^X 



(6) 



tanh[£ k _ q/2 /(2fc B T)] + tanh[£ u+q/2 /(2fc B T)] 



Ck-q/2 

so that one obtains 
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S b (q,0) = E,(0,0) 
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for small q with any anisotropy factor compatible with 
the point-group symmetry of the cuprates. Here M* is 
the boson mass, calculated analytically in Ref.^1 with 
the isotropic exchange interaction and parabolic fermion 
band dispersion (see also Ref.0), and h = 1. The BCS- 
like equation (2) has a nontrivial solution for A k at T = 
T r , if 



Eo = -£6(0,0). 



(8) 



Substituting Eq.(7) and Eq.(8) into the sum-rule, Eq.(5) 
one obtains a logarithmically divergent integral with re- 
spect to q, and 



„ const 

T R = -^—^ = 0. 
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The devastating result, Eq.(9) is a direct consequence 
of the well-known theorem, which states that BEC is im- 
possible in 2D. This is true for non-interacting bosons. 
Remarkably it appears also true for bosons hybridized 
with fermions in 2D, Eq.(9). Any dynamic repulsion be- 
tween bosons could provide an infrared cut-off of the in- 
tegral in Eq.(9) resulting in a finite critical temperature 
in 2D. But it has to be treated beyond MFA level USUI, 
and in no way does it make MFA El 

a meaningful 

approximation. 



One may erroneously believe that MFA results can 
be still applied in three-dimensions. However, increas- 
ing dimensionality does not make MFA a meaningful ap- 
proximation either. This approximation leads to a naive 
conclusion that a BCS-like superconducting state oc- 
curs below the critical temperature T c ~ /iexp (—Eo/z c ) 
via fermion pairs being virtually excited into unoccupied 
virtual bosonic states [9J, |10J, [l2| . Here z c = g 2 N(0) and 
N (0) is the density of states (DOS) in the fermionic band 
near the Fermi level \x. However, the Cooper pairing of 
fermions is not possible via virtual unoccupied bosonic 
states in 3D BFM either. Indeed, Eq.(8) does not de- 
pend on the dimensionality, so that the analytical con- 
tinuation of Eq. (4) to real frequencies w yields the par- 
tial boson DOS as p(w) = (1 - 2n B )S(uj) at T = T c for 
q = in any-dimensional BFM and for any coupling with 
fermions. Hence, the Cooper pairing may occur only si- 
multaneously with the Bose-Einstein condensation of real 
bosons in 3D BFM 13] . The origin of the simultaneous 
condensation of the fermionic and bosonic fields in 3D 
BFM lies in the complete softening of the boson mode at 
T = T c caused by its hybridization with fermions. 

Taking into account the boson damping and dispersion 
shows that the boson spectrum significantly changes for 
all momenta. Continuing the self-energy, Eq.(6) to real 
frequencies yields the damping (i.e. the imaginary part 
of the self-energy) as |l3j 



7(q,w) 



In 



cosh(<z£ + u;/(4fc B T c )) 



cosh(-g£ + w/(4fc B T c )) 



(10) 



where £ = Ui?/(4fc B T c ) is a coherence length. The damp- 
ing is significant when q£ « 1. In this region 7(q, ui) = 
u)irz c /(8kBT c ) is comparable or even larger than the bo- 
son energy u. Hence bosons look like overdamped diffu- 
sive modes, rather than quasiparticles in the long-wave 
limit |j3l Ho ], contrary to the erroneous conclusion of 
Ref. [ll| , that there is 'the onset of coherent free-particle- 
like motion of the bosons' in this limit. Only outside the 
long-wave region, the damping becomes small. Indeed, 
using Eq.(10) one obtains j(q,co) — u>irz c /(2qvF ) << u, 
so that bosons at q >> l/£ are well defined quasiparti- 
cles with a logarithmic dispersion, uj(q) = z c ln(q^) [l3| . 
As a result the bosons disperse over the whole energy 
interval from zero up to Eq. 

The main mathematical problem with MFA stems from 
the density sum rule, Eq.(5) which determines the chem- 
ical potential of the system and consequently the bare 
boson energy Eo(T) as a function of temperature. In 
the framework of MFA one takes the bare boson en- 
ergy in Eq.(2) as a temperature independent parame- 
ter, Eo = z c ln(fi/T c ) 12], or determines it from the 
conservation of the total number of particles, Eq.( 3) ne- 
glecting the boson self-energy E b (q,f2) 0, E3 [3 111). 
Then Eq.(2) looks like a conventional Ginzburg-Landau- 
Gor'kov equation |l9j linearized near the transition with 
a negative coefficient a cx T — T c < at T < T c , 



aA(r) = 0, 



(11) 
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where 



a = 1 



£ b (0,0) 

Eq 
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|tf(0) m 
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Hence, one concludes that the phase transition is almost 
a conventional BCS-like transition, at least at Eq 3> T c 
|fll lT(ilT^ . Also, using the Gor'kov expansion in powers 
of A in the external field, one finds a finite upper critical 
field H c2 (T) [H|. 

However, these findings are mathematically and phys- 
ically flawed. Indeed, the term of the sum in Eq.(5) with 
fi„ = is given by the integral 



T 



27r3^ + E b (q,0)' 



(13) 



The integral converges, if and only if Eq ^ — Eb(0,0). 
Hence the coefficient a(T) in Eq.(ll) can not be negative 
at any temperature below T c contrary to the MFA result 
jl2l Il5| , which violates the density sum-rule predicting a 
wrong negative ot(T). 

Since a(T) ^ 0, the phase transition is never a BCS- 
like second-order phase transition even at large Eq and 
small g. In fact, the transition is driven by the Bose- 
Einstein condensation of real bosons with q = 0, which 
occur due to the complete softening of their spectrum 
at T c in 3D BFM. Remarkably, the conventional upper 
critical field, determined as the field, where a non-trivial 
solution of the linearised Gor'kov equation 0] in the 
external field occurs, is zero in 3D BFM, H C 2(T) = 0, 
because a{T) ^ below T c (for more details see [2ljV 
It is not a finite H c2 (T) found in Ref. [lj| using MFA. 
Of course, like in 2D BFM, the dynamic repulsion be- 
tween bosons could provide a finite H C %{T) in 3D BFM, 
similar to the case of intrinsically mobile 3D bosons [2(| . 
However, to provide a finite H C 2(T) the dynamic repul- 
sion has to be treated beyond MFA, as discussed in Ref. 
[2"o|. Even at temperatures well below T c the condensed 
state is fundamentally different from the BCS-like MFA 
ground state, because of the pairing of bosons |2l|. It 
is similar to a pairing of supracondensate helium atoms 



in 4 iJe, proposed as an explanation for a low density 
of the single-particle condensate^. The boson pairing 
appears in 3D BFM due to the hybridization of bosons 
with fermionic condensate. It is not expected in the 
framework of MFA, where the effective interaction be- 
tween bosons was found repulsive by integrating out the 
fermionic degrees of freedom with an incorrect (constant) 
Eq The pair-boson condensate should significantly 

modify the thermodynamic properties of the condensed 
BFM compared with the MFA predictions. 

This qualitative failure of MFA might be rather un- 
expected, if one believes that bosons in Eq.(l) play the 
same role as phonons in the BCS superconductor. This 
is not the case for two reasons. The first one is the den- 
sity sum- rule, Eq.(5), for bosons which is not applied to 
phonons. The second being that the boson self-energy 
is given by the divergent (at T = 0) Cooperon diagram, 
while the self-energy of phonons is finite at small cou- 
pling. 

I conclude that MFA results for the boson-fermion 
model do make any sense neither in two nor in three 
dimensions, because the divergent self-energy has been 
neglected in calculating T c and H C 2(T). Any repulsion 
between bosons, induced by hard-core effects and/or by 
hybridization with fermions could not make MFA results 
U IH IH 111 El meaningful. It is well known 
that the " infrared-save" 2D theory by Popov and oth- 
ers 0, is actually not a mean-field theory. In 3D 
BFM the phase transition is due to the Bose conden- 
sation of real bosons rather than virtual ones. Hence, 
the mean-field theory cannot be applied for a description 
of 3D BFM either. There is no BCS-like state in any- 
dimensional BFM, and no BCS to local pair crossover. 
The common wisdom that at weak coupling the boson- 
fermion model is adequately described by the BCS the- 
ory, is negated by our results. 
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